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DEFLECTION AND STRESS ANALYSIS OF THIN SOLID WINGS OF ARBITRARY PLAN FORM 

WITH PARTICULAR REFERENCE TO DELTA WINGS 1 

By Manuel Stein, J. Edward Anderson, and John M. Hedoepeth 


SUMMARY 

The structural analysis oj arbitrary solid cantilever wings by 
small-deflection thin-plate theory is reduced to the solution of 
linear ordinary differential equations by the assumption that 
the chordwise deflections at any spanwise station may be ex- 
pressed in the form of a power series in which the coefficients are 
functions of the spanwise coordinate. If the series is limited to 
the first two and three terms (i that is, if linear and parabolic 
chordwise deflections, respectively, are assumed), the differential 
equations for the coefficients are solved exactly for uniformly 
loaded solid delta wings of constant thickness and of symmetrical 
double-wedge airfoil section with constant thickness ratio. For 
cases for which exact solutions to the differential equations 
cannot be obtained, a numerical procedure is derived. Experi- 
mental deflection and stress data for constant-thickness delta- 
plate specimens of 45° and 60° sweep are presented and are 
found to compare favorably with the present theory. 

INTRODUCTION 

One of the present trends in the development of high-speed 
airplanes and missiles is toward the use of thin Iow-aspect- 
ratio wings. The structural analysis of these wings often 
cannot be based on beam theory since the structural defor- 
mations may vary considerably from those of a beam and, 
indeed, may more closely approach those of a plate. In 
cases where the wing construction is solid or nearly solid the 
use of plate theory in the analysis is particularly valid, and 
it is this type of wing which is considered in the present 
report. 

Exact solutions to the partial-differential equation of plate 
theory are not readily obtained, especially for plates of 
arbitrary shape and loading; however, a number of approxi- 
mate solutions to specific problems on cantilever plates have 
appeared in the literature (see, for example, refs. 1 to 7). 
Of the approaches used in these references, only the one in 
references 6 and 7 is readily applicable to plates of arbitrary 
plan form, thickness distribution, and load distribution; thus 
it is the most useful one for the analysis of actual wings. 

In reference 6 the cantilever-plate problem is simplified by 
the assumption that the deformations of the plate in the 
chordwise direction (parallel to the root) are linear. By 
minimizing the potential energy of the plate, the partial- 
differential equation of plate theory is replaced by two 


simultaneous ordinary differential equations for the spanwise 
variations of the bending deflection and twist. In reference 7 
the same ordinary differential equations are obtained in a 
different manner. Refinement of the analysis by inclusion 
of the effect of parabolic, cubic, or higher-order chordwise 
camber terms is indicated in reference 6, and as the order of 
refinement is increased a corresponding increase in the num- 
ber of ordinary differential equations is obtained. 

In the present report, which is an extension of reference 6, 
a general set of ordinary differential equations is presented 
which may be used to obtain any desired degree of approxi- 
mation to the deflection of the plate. These equations are 
solved exactly for several cases of delta plates under uniform 
load first by considering linear chordwise deformation only 
and second by including the effect of parabolic chordwise 
camber. Comparisons are drawn between the stresses and 
deflections computed from the equations of each approxi- 
mation and also with some experimental results. 

The differential equations presented contain coefficients 
that depend on the plan form and stiffness distribution of the 
plate and on the loading. In this report, the plates con- 
sidered in detail have coefficients such that the differential 
equations can be solved exactly; however, in cases for which 
exact solutions cannot be obtained a numerical procedure 
must be used. One such procedure is derived and its 
accuracy is demonstrated. 


I 

c 

V 

t 

tao 

D 

T) 


E 


v- 

w 


y> * 


SYMBOLS 


length of plate measured perpendicular to root 
root chord of plate 

lateral load per unit area, positive in in- 
direction 

local thickness of plate 
average thickness of plate 


local flexural stiffness, 


Et 3 

12(1 m 2 ) 


flexural stiffness based on average thickness 
Et 3 


12(1 M 2 ) 

modulus of elasticity of material 
Poisson’s ratio 


deflection of plate, positive in ^-direction 
coordinates defined in figure 1 


1 Supersedes NACA TN 2621, “Deflection and Stress Analysis of Thin Solid Wings of Arbitrary Plan Form With 

Anderson, and John M. Hedgepeth, 1952. 


Particular Reference to Delta Wings” by Manuel Stein, J. Edward 
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Figure 1. — Coordinate system used in the present analysis for a canti- 
lever plate of arbitrary shape with arbitrary thickness variation. 


<Pn 

function of x, coefficient in power series for 


N 

deflection w='^ l <p n (x)y n 

n =0 

ci(x), c 2 (x) 

functions defining plan form (see fig. 1) 

Xi 

X 

variable obtained by transformation Xi = \ — j 

<T X ) Gy 

normal stresses 

Txy 

shear stress 

G 

maximum principal stress 

X 

aspect-ratio parameter, (1 ~~m) 


RESULTS 

The derivation of the general set of ordinary differential 
equations is given in appendix A. The general procedure 
outlined in reference 6 is followed; that is, the deflection 
of the plate w is expanded into a power series in y the chord- 
wise coordinate with coefficients which are functions of x 
the span wise coordinate (see fig. 1) 

w=(p 0 (x) + ipi(x)y + ip2(x)y 2j r • • • + ( 1 ) 

Equation (1) is substituted into the expression for the po- 
tential energy of the plate-load combination which is in 
turn minimized by the calculus of variations with respect 
to each of the coefficients <p n . The results of the variational 
procedure appear as JV+1 simultaneous differential equa- 
tions with the coefficients <p n as unknowns. 

By taking a sufficient number of terms in tin* expansion 
of w, the resulting differential equations can be used to ob- 
tain any desired degree of accuracy in the solution for the 
deflections of any given cantilever plate subjected to an 
arbitrary lateral load. Of most interest, perhaps, are the 
particular cases for N= 1 and N=2, which are obtained 
from the general set of equations and are simplified in 
appendix A. The case for N= 1 (also derived in refs. 6 


and 7) includes linear chordwise deflections, and the case for 
N=2 takes into account parabolic chordwise curvature. 
Although for most practical problems the solutioa by the 
parabolic theory should be adequate, cases might exist in 
which cubic, quartic, or even higher-order chordwise terms 
should be included, depending on the convergence of the 
series for the configuration considered. 

The particular equations for N=l and N= 2 are used to 
determine the deflections and stresses of the following can- 
tilever plates subjected to uniform lateral load: 

(1) A 45° delta plate of uniform thickness 

(2) A 60° delta plate of uniform thickness 

(3) A 45° delta plate of symmetrical double-wecge airfoil 
section with constant thickness ratio 

Fortunately, for these configurations, the solution can be 
carried out exactly by both the linear and parabolic theories, 
and the details of these exact solutions are included in 
appendix B. In general, however, exact solutions cannot 
be obtained and some numerical method must be used. 
One such method, based on replacing derivatives by their 
first-order-approximation difference forms, is derived in 
appendix C. 

A summary of the results for the three particular problems 
is shown in figures 2 to 11. Deflections obtained by the 



Figure 2. — Deflections of a 45° delta plate of uniform thickness 

uniform load. 
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linear theory and the parabolic theory for the three con- 
figurations are compared in figures 2, 3, and 4. Stresses 
obtained by the linear theory and the parabolic theory 
for the three configurations are compared in figures 5, 6, 
and 7. Where available, experimental deflections and 
stresses are also shown in these figures. The details of the 
procedure used to obtain the experimental deflections of 
the 45° and 60° uniform-thickness plates and the experi- 
mental stresses in the 45° uniform-thickness plate are con- 
tained in appendix D ; whereas the experimental root stresses 
for the 60° uniform-thickness plate were obtained from 
reference 8. Figures 8 to 11 present the comparison be- 
tween deflections and stresses computed from the exact 
solutions of the differential equations and those computed 
from the numerical solutions of the same equations. 

DISCUSSION 

The results shown in figures 2 and 3 indicate that, with 
regard to deflections, either the linear theory or the parabolic 
theory is adequate for the case of a constant-thickness delta 
plate subjected to a uniform load, the comparison being some- 
what better for the 60° plate than for the 45° plate. If 
accurate slopes in the chordwise direction (angle of attack) 



Figure 3. — Deflections of a 00° delta plate of uniform thickness under 

uniform load. 


are desired, however, the parabolic theory must be used 
because the error in the angle of attack as computed by the 
linear theory is as much as 30 percent (see figs. 2 and 3). 
The appreciable anticlastic curvature, evidenced by the 
experimental results of figures 2 and 3, may be important 
aerodynamicallv and is, of course, not taken into account by 
the linear theory. 

The apparent convergence of the aforementioned series 
in the case of the double-wedge-section plate (see fig. 4) 
implies that the linear theory is adequate for this case. The 
lack of chordwise curvature in the result obtained by the 
parabolic theory is attributable to the fact that the natural 
tendency of the plate to have anticlastic curvature is canceled 
by the opposite tendency of the thin edges to bend down 
under the load. Unfortunately, no experimental results 
are available for this configuration. 

In figure 4 the plate stiffness D in the nondimensional 
parameter wD/pl 4 is the local value of D at a point where 
the thickness is equal to the average thickness of the plate 
as a whole. Thus the results of figure 4 are comparable 
with the results of figure 2 on an equal-weight basis. It 
can be seen that the deflections of the double-wedge-section, 
constant-thickness-ratio plate are everywhere less than 



Figure 4. — Deflections of a 45° delta plate of symmetrical double- 
wedge airfoil section and constant thickness ratio under uniform load. 
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those of the uniform-thickness plate although they increase 
rapidly near the tip. This curling-up or singularity in slope 
at the tip is a result of using a small-deflection theory and 
probahlv would not be so marked in an actual case. 

The stress results for the 45° and 60° uniform-thickness 
delta plates indicate that both the linear and the parabolic 
theories are adequate and that the parabolic theory is better 
than the linear theory only near the root. It should be noted 
that, although the maximum principal stress over a large 
part of the 45° plate is plotted in figure 5, only the stresses 

normal to the root along the line j= 0.0087 of the 60° plate 

are plotted in figure 6 since only these stresses are given in 
reference 8. (The maximum principal stress and the stress 
normal to the root are theoretically equal at the root since 
the root shear stress is zero.) 

Experimental data are lacking for the double-wedge-section 
delta plate and, therefore, only theoretical stresses are shown 
in figure 7. As in the case of deflections, the results obtained 
from the linear theory and those obtained from the parabolic 
theory are almost coincident, the difference being greatest 
near the root. Figure 7 has also been plotted so that the 
results are directly comparable with those for the 45° 



Figure 5. — Maximum principal stresses in a 45° delta plate of uniform 
thickness under uniform load. 


uniform-thickness plate in figure 5 on an equ il-weight 
basis. As can be expected, the double-wedge-section, 
constant- thickness-ratio plate is a better design stn cturally; 
the stresses in the double-wedge-section plate are everywhere 
smaller and are almost constant in the span wise direction. 

The theoretical results in figures 2 to 7 have been obtained 
from exact solutions of the differential equations of the 
linear and parabolic theories. In order to test the reliability^ 
of the numerical method derived in appendix C, the differ- 
ential equations were also solved numerically. The results 
shown in figures 8 and 9 indicate that the agreement is good 
between the numerical solution in which five equal intervals 
were used and the exact solution of the differential equations 
for the case of the 45° uniform-thickness plate. The same 
good agreement can be expected in other cases where the 
thickness and load distributions are not too erratic and where 
the plate stiffness does not go to zero at the tip — that is, 
when no singularities appear at the tip. 



Figure 6. — Normal-stress distribution near the root ^at -=0.0087^ of 
a (>0° delta plate of uniform thickness under uniform load. 
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Since the efficacy of the numerical method depends on 
how well parabolic arcs fit the various functions between 
stations, serious error can result from blind application. An 
example of the seriousness of these errors and of the manner 
in which they can be remedied is shown in figures 10 and 11. 
In these figures a comparison is made between exact and 
numerical results obtained on the 45° double-wedge-section, 
constant- thickness-ratio plate. As can be expected, the 
five-station numerical solution fails to follow the exact solu- 
tion in the neighborhood of the singularity at the tip. Since 
the region of trouble is localized at the tip, a reasonable 
remedy would be to decrease the spacing of the station points 
near the tip. This decrease in spacing may be accomplished 
either by using a greater number of equally spaced stations 
or by using unequally spaced stations crowded near the tip. 
The increase in accuracy obtained by increasing the number 
of equally spaced station points to ten is shown in figures 
10 and 11. 

CONCLUDING REMARKS 

The general method presented herein for finding deflec- 
tions and stresses of solid or nearly solid wings is, in principle, 



Figure 7. — Maximum principal stress in a 45° delta plate of symmet- 
rical double-wedge airfoil section and constant thickness ratio under 
uniform load. 


capable of yielding arbitrarily accurate results for any con- 
figuration. It is seen that, for the examples considered, 
only the first two or three terms in the series expansion need 
be considered to obtain adequate accuracy. In addition, 
for most practical plate-like wings with clamped roots the 
first two or three terms will probably be adequate, although 
problems may exist wherein more terms are needed. 

The numerical procedure, derived for application in cases 
where exact solutions cannot be obtained, gives good agree- 
ment when compared with exact solutions if enough stations 
are taken along the span. The necessary number of stations 
is dependent on the type of thickness and loading distribution 
considered, five equally spaced stations being enough for 
the uniform-thickness delta wing subjected to uniform 
loading and ten being necessary for the double-wedge-section, 
constant- thickness-ratio delta wing subjected to uniform 
loading. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November SO , 1951 . 



Figure 8. — Numerical and exact solutions of the differential equations 
for the deflections of the free edges of a 45° delta plate of uniform 
thickness under'uniform load. 
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Figure 9. — Numerical and exact solutions of the differential equations 
(obtained by assuming linear chordwise deflections) for the maximum 
principal stresses along the free edges of a 45° delta plate of uniform 
thickness under uniform load. 





Figure 10. — Numerical and exact solutions of the differential equations 
(obtained by assuming linear chordwise deflections) for he deflec- 
tions along the free edges of a 45° delta plate of symmetrical double- 
wedge airfoil section and constant thickness ratio under uniform load 
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Figure 11. — Numerical and exact solutions of the differential equa- 
tions (obtained by assuming linear chord wise deflections) for the 

maximum principal stress along the line of a 45° delta 

plate of symmetrical double-wedge airfoil section and constant 
thickness ratio under uniform load. 


APPENDIX A 


DERIVATION OF DIFFERENTIAL EQUATIONS 


The structure considered herein is a thin, elastic, isotropic, 
cantilever plate of arbitrary plan form and slowly varying 
thickness subjected to distributed lateral load (see fig. 1). 
By assuming that the deflection of the plate can be repre- 
sented by a power series in the chord wise coordinate and 
by applying the minimum-potential-energy principle, a set 
of ordinary differential equations in the spanwise coordinate 
is obtained from which the coefficients of the power series 
may be determined. First the general set of equations is 
derived; then the particular equations for the cases of linear 
chordwise deflections and parabolic chordwise deflections are 
deduced from the general set and simplified. 


General equations. — The potential energy of the systeir 
under consideration is 


Potential ener 


-n 


c - (r) ( D(x,y) /b 2 w\ 
•>» ( 2 ' 1 1 ~ 2 J + 


/ b 2 w V 

wv 


+ 


b 2 w b 2 iv 


bx 2 by 


3 + 2 ( 1 -*) 


LV C )X 2 J 

( bx by ) ]~P( X >^ W ) d V dx ( Al ) 


in which 


D(x,y) 


E[ K^,y)] 3 

12(1 -m 2 ) 


and p(x,y) is the distributed lateral load. 


266368—54 2 
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The assumption is made that the deflection w can be represented by tbc power series 

N 

W=S <Pn(x)y n 
77=0 

Substitution of this expression for w into equation (Al) gives 

ri (i N N 

Potential energy = dxl-^z X 2 [«m+«+i^m / Vn // + M (^~1) (n— 1) a m+n _ 3 ^ m (p n + 

J0 ( ^ 771=0 71=0 

2 1X71 (?l + \<Pm ¥ ? n - i - 2( 1 /x)/W-7ld m + n _ V 3 // ] + 

n=o ; 

in which 

<*,= f D(x,y)y r ~ l dy (r= 1 , 2 , .. . 2.Y+ 1 ) ' 

J Cl(x) 

J * c 2 (a:) 

, p(x,y)y T ~ l dy (r= 1,2, . . . V+l) 

and the primes denote differentiation with respect to x. 

Minimization of the potential energy by means of the calculus of variations gives 

^(Potential energy) =0 


(A 2) 


(A3) 


(A4) 


ri ( j N N 

= I X* X* \Q'm + n + \(}Pm <Pn d<Pm )“f "77171(771 1) ( 71 \)(l m + n -z(ip m 5(p n T <PnS*Pm) T" 

Jo ( - 771=0 71 = 0 

2 fi7i(n— l)a m+n _i(^ m "6^ n + ^ w 6^ m ") + 2(l — /x)mna m+n -.i(^ m , 5v? n ' + ^> w '5^ m / )]“S P*+ 1 \ 

71 = 0 ) 

Integrating by parts and collecting terms results in 

0= ( (Ix^Sipn js [(a m+ #+lV’n.")" + M7»(»7— 1) («U + n- — 2 (1 — /<) »i» (am + B-UPmO' + /*«■ (»~ 1) am+n-HPm" + 

Jo 77=0 ( 777=0 

W A r V y l N N 

yZ^iPn X [^m + n + 1^7«" + M m ~ 1) O-m + n-l^m] ? ) X XI [(«//! + ti + 1^/7, 'O' + 

77 = 0 7/7 =0 ) 0 l 77 = 0 777=0 

fX Vfl (jYl 1 ) (tt m -f n — 1 *P /«) ^(l /x) W /7 d m + /i - 1 ^ /// ] | ^ ^ 

Everywhere in the region of the plate, except at the boundary j= 0, the variation of w is arbitrary. At z = 0 the 
cantilever boundary conditions 


yield 

*.(0) = *.'(0)=0 ( 71 = 0 , 1, . . . N) (A6) 

and therefore the variation in these quantities must also be zero at x=(). 

Equation (A5) is then satisfied if, in addition to equation (A6), 

X [(a m + n + 1 V?m")" + M^(^— l)(a w + n _i^ m )" — 2(1— rimnCdn + n-ipmO' + MnCn— 1 )a m+n - 1 <p m " + 

777 = 0 

mn(m—\)(n— l)a m+n _ 3 (f> m \ = p n+ i (n= 0 , 1 ,... A) (A7) 

ii2[a m+n+ i<p m ''+iim(7ri-l)a m+n -i<p m ]x^i=0 (n= 0, 1, . . . N) (A8) 

7/7 = 0 

and 

X) [( a rn + n + l < ^/« // ) /- l _ A tm ( m — 1) (Pm + n- \<P hi)' 2(1 ix)m7l(l m + n -\(f> m ] x =l = 0 (ft=0, 1, ... A) (A9) 

m=0 

Equations (A7) form a set of N+\ simultaneous ordinary differential equations for the functions ^„(x). The fund ions <p n 
are completely determined by these diff erential equations and the boundary conditions (A6), (A8), and (A9). 
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Particular case of linear chordwise deflections. — If 

A r =l, the deflection function becomes 

w — < Po J ry ( P\ (A 10) 

a linear function in the chordwise direction, where ? 0 is the 
bending deflection and <p\ is the twist. Equations (A7) 
become 

(<W ')"+ (a 2 <Pi")" = P 1 (Al 1) 

((i 2 <Po")" -h — 2(1 — aO ( Q’\<Pi)' = p 2 (A 12) 

The root boundary conditions, given by equation (A6), 
become 

<A) (0) = w (0) = <p\ (0) = ipi ' (0) = 0 (A 1 3) 

The tip boundary conditions, given by equations (A8) and 
(A9), become 

(®lW / +^2^1 /, )ar-I=0 (A 14) 

( a 2<Po' ' + Q'3 <Pl ' ') x= l = 0 (A 1 5) 

[{am , y+(a 2 i Pl ' , y] x=l ^{) (AlG) 

[{<h<Po")' — v)<Li<pi'] x =i=0 (Al7) 

Equations (All) to (A17) arc the differential equations and 
corresponding boundaiy conditions presented in reference 
0 (if only distributed load is considered) where the symbols 
IT and 0 are used instead of ? 0 and <p u respectively. 

If equation (All) is integrated twice and the boundary 
conditions (A14) and (AlG) are used, 

^°"=— ~ <Pi"+-~ L f \ pi dx 2 (A 18) 

a \JrJz 

Substitution of <p 0 " into equations (A 12), (A 15), and (A 17) 
gives 

2(1 — n) (ai<p/)'=p 2 — (— J* J p t dx 2 ) (Al9) 

(bm")x=i=0 (A20) 

[{b lVl "Y -2(1- ^a m '] x=l =0 (A21) 

in which 

? a 2 2 

b\ — a s — 

d\ 

If equation (A 19) is integrated once and the boundary 
condition (A21) is used, 

(bnpi')'~ 2(1— ii)anpi= — J p 2 rfz — J Pi dx 2 ^ (A22) 

The differential equation (A22) is a second-order differential 
equation in ?/. The twist ?i and then the bending deflection 
<p 0 are obtained by solving equations (A22) and (Al8), respec- 
tively, by applying the boundary conditions (Al3) and (A20). 

Particular case of parabolic chordwise deflections. — The 
effect of parabolic chordwise camber may be included by 


letting N= 2 in the general power series (eq. (A2)). If N=2, 
the deflection function becomes 

W— <A) + 2/<Pl + 2/V2 

Here ? 2 represents the spanwise distribution of parabolic 
chordwise camber. For this case the differential equations 
(A7) become 

(<WT + (&2<pi // ) // 4 - (&z(p 2 f )” -\-2n(d\<p^)' f ~P\ (A23) 

) " T (dz<P\ ") " + (tf'4^2 ") " T 2 M (d 2 <p 2 ) " — 

2 (l-n)[(a m 'y +2 (a 2 <p 2 ) ' J =p 2 ( A24 ) 

(a 3 ? 0 ")" ~\r (d 4 ip\ ")*-{- (d 5 (p 2 " ) " + 2 n [di (p 0 " + d 2 ip x " + a 3 ? 2 " + 

( 0 ^ 2 ) "] 4 (1 — fj.) [(d 2 <p x ') ' +2 (d 3 <p 2 ) '] +4ai? 2 ~Pz ( A25) 

with the boundary conditions 

^0 (0) = ? 0 ' (0) = ?! (0) = ?! ' (0) = ?2 (0) - tpt’ (0) = 0 ( A2G) 

( a \W' -\~(i2<p 1 ' +a 3 ? 2 ,r T2/zai? 2 )x=z — 0 (A27) 

(o, 2 ipo" -j- a 3 ?i n T d 4 (p 2 " T 2jjLd 2 <p 2 ) 3 = j = 0 (A28) 

+ &4<Pl ” + dbVl' + 2 jLia 3 ? 2 ) x= / = 0 (A29) 

[(ai? 0 ")'+ (a 2 ?i // )'+ (a 3 ? 2 ,, ) , + 2 ju(a 1 ? 2 )']z =/=0 (A30) 

[(a 2 ? 0 /, ) / + («4^2 /, ) / + 2/i(a 2 ?2) / -“ 

2 (1 — m) ( a i <P\ ' + 2a 2 ? 2 / )]x= 1 ~ 9 (A3 1 ) 

[(&3^</ , ) / + (tt 4 ?l /, ) , + (0-5tf , 2 // ) , + 2p(a 3 ? 2 ) / — ' 

4(1 (^2^i / + 2a 3 ? 2 / )]x=/— 0 (A32) 

If equation (A23) is integrated twice and the boundary 
conditions (A27) and (A30) are used, 

^o"=— ^ ?/'— “ ? 2 " — 2 /x?2+— f ( Pi dx 2 (A3 3) 

®l J.r Jx 

Substitution of ? 0 " into the remaining differential equations 
and boundary conditions results in 

(WT +(W')''-2( 1 -/x) [(a 1 ? 1 ') , + 2 (a 2 ? 2 ')'] 

&, y <A34) 

(6 2 <^i ) ~\~(b 3 ip2 )" 4(1 m) [(^2^i / ) / T2(a 3 ? 2 / ) / ]T4(l — p?)d\<p 2 

*7' (A35) 


(&i^i ,, 4"62^2 ,, )*-i == 0 (A3 6) 

[(6i?i , 0 , + (W0 , -2(l~M)(ai? 1 ' + 2a 2 ?2 , )]x»*=0 (A37) 

(62?/' + b 3 <p 2 ") x =i=0 (A38) 

[(62^1 ) T(^3<^2 / ) / 4(1 m) (d2(p\ -b2a 3 ? 2 , )]x=j == 0 (A39) 

?1 (0) = ^ ' (0) = ? 2 (0) — ? 2 '(0)=0 (A40) 
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in which 


i a 2 

b ' =<h ~, h 


b 2 =(i\— 


a 2 (h 


6 3 =« 5 - ( - 


If equation (A34) is integrated and the boundary condition 
(A3 7) is used, 


Stresses. -After the approximate deflection of the plate 
is determined from equations (A 18) and (A22) or from equa- 
tions (A 33), (A3 5), and (A41), the extreme-fiber stresses may 
be calculated from the well-known equations of tl in-plate 
theory, which are (see, for example, ref. 9): 

67 ) (b 2 w b 2 w\ 

Wv 

6D /d 2 w d 2 w\ 
a “~ t * \dy 2+ ^ dY 2 ) 


{b (b 2 <pi"Y — 2(1 — n) -\-2a 2 >pY) 

(A41) 

Thus <p x and cp 2 are obtained by solving equations (A35) and 
(A41) with the boundary conditions (A36), (A38), (A39), and 
(A40) . Subsequently, <? 0 can be obtained by solving equa- 
tion (A33) with the boundary conditions <p 0 (0) = <Po'(0) =0* 


6 ( 1 — n)D b 2 w 
Txv ~~ t 2 bx by 

The maximum principal stress cr at any point in the plate 
can be determined from 


O'* 


+ 0V 


± 2 V(°*x“ 


-<r„) 2 +W 


2 


APPENDIX B 

EXACT SOLUTIONS OF DIFFERENTIAL EQUATIONS FOR SOME SPECIFIC DELTA-PLATE PROBLEMS 

chordwise deflections result : 


The differential equations of appendix A for linear and 
parabolic chordwise deflections are solved exactly for uni- 
formly loaded delta plates of constant thickness and of sym- 
metrical double-wedge airfoil section with constant thickness 
ratio. The equations for deflections obtained by the linear 
theory are presented in terms of the aspect-ratio parameter 
X for both kinds of delta plates. The equations for deflec- 
tions obtained W the parabolic theory are presented for 

l /3 i 

-=1 and with n=- for the constant-thickness delta plate 

Co 6 

and for - = 1, also with n=h> for the delta plate of symmet- 

C o 

rical double-wedge airfoil section with constant thickness 
ratio. 

If the z-axis is passed through the edge perpendicular to 

x 

the root and the substitution Xi=\—j is made, the differen- 
tial equations are clearly of the homogeneous type for which 
the solutions are of the form Xi y , where 7 is a constant. For 
the configurations considered, the functions that define the 
plan form (see fig. 1) are then Ci(z) = 0 and c 2 (x) = cx \ , where c 
is the root chord. In all the equations of this appendix the 
primes denote differentiation with respect to the new inde- 
pendent variable X\. 

DELTA PLATE OF UNIFORM THICKNESS UNDER UNIFORM LOAD 

Since the stiffness D is a constant for uniform-thickness 


(A4)) become 


in 

the differential 


_Dc n 

X\ n 

a n 



n 


a 2 2 _ 

Dc 3 , 

= a 3 

a>\ ~ 

= _ 12~ X ' Z 

• n 

d 2 d% 

Dc* 

— / y* - ** 


d\ 

12 

n . - 

_al_ 

4 Dc 5 . 

7* , 0 

a o 

d\ 

45 1 


n - V —x « 
Vn~ n 


(Bla) 

(Bib) 

(Blc) 

(BUI) 

(Ble) 


Solution for linear chordwise deflections. — If the co- 
efficients given by equations (Bl) are substituted into equa- 
tions (A22) and (Al8) and the independent variable is 


changed to Xi=l —j-> the following 


equations for linear 


where 


(xi 3 (pi")' — 1 6\ 2 x t <f!i' = — 2 x , 3 


„ C „ . p / 4 2 

<Po — — 2 Xl(pl QD Xl 


(B2) 

(B3) 


X=- 


V 2 


(l— M) 


The boundary conditions to be used with these equations are 
obtained from equations (A 13) and (A20) and are 


^o(l) = ^o'U)=^i(l) = ^i'(l) = 0 
(Xi 3 <P\")x 1= o = 0 

The general solution of equation (B2) is 


where 




(B4) 

(B5) 

(BO) 


7=^1 + 16X 2 


and Ai and A 2 are arbitrary constants. Since X 2 is inherently 
positive, the boundary condition (B5) requires that ^1 2 =0. 
One integration of equation (B6) and the application of the 
conditions <pi(l) = <p/(l) = 0 yields 


<P 1 = 


1 pP / X 1 7 — 1 _Xi 3 

~4(i — 2\ 2 ) TTc : 


') 


(B7) 


If equation (B3) is solved for <p 0 with the conditions 
V?0 (1) = ^ 0 '(1) = 0, the result is 

"-fjj r=W[f (5 - 4X ’) O—- ‘ t 1 *)- 


— ( 1 


0 


-Xi — 


1 — Xi y+l 

7 + 


f)] 


(B8) 


Substitution of equations (B7) and (B8) into the equation 

gives the expression for the deflection w of the plate under 
the assumption of linear chordwise deflections. 

Solution for parabolic chordwise deflections. — If the co- 
efficients given by equations (Bl) are substituted into equa- 
tions (A41), (A35), and (A33) and the independent variable 

x 

is again changed to Xi=l—-jy the following equations for 


11 


12 


REPORT 1131 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


parabolic chord wise deflections result: 


Zl 


(B9) 


(x.V 1 T , + j5(* I 5 ^2 , 0 ,, -16X 2 [(x 1 V. , )'+|(^. 3 ^2 , )']+ 

64 . 1+M 4 / 2pX 2 \ pf 4 , 

t x i-M WJ= -H' + r^)^ 3 '' (B10) 

V'=-| x,V'-2 M /V 2 +^^ (B 11) 

The boundary conditions to be used with these equations are 
*(l) = *'(l) = Ml) = w'(l) = *»(l) = «'(l)=0 (B12) 

faV'+ahWOft-o^O (B 1 3) 

(x 1 V 1 "+j|^ , W , )_ o =° (Bl 4) 

[(siVO'+jf (mVO'-iex* (x.V.'+l 

(Bl 5) 

The homogeneous solutions of the simultaneous equations 
(B9) and (BIO) are of the form 

<p x '=Axi y ~ l 
(P‘2 = J^X\ y 1 

Substitution of these expressions into the homogeneous parts 
of equations (B9) and (BIO) leads to the following charac- 
teristic equation from which X may be determined: 

7 6 — 6(1 + 1 6X 2 )7 4 +[^320 ( 4 + J ^ ) X 4 + 480X 2 + 9 Jy 2 — 

4 [ i28 ° x"+ 8 °( 4 +|:;) x 4 +96x 2 +ij=o 

(B 1 6) 

and gives the following relationship between A and li : 

4 ( iT(t- 2)^+1 )- lfi X 2 l,R 

-4 = -( 7 - 1 ) L 7 2“ fZ+ 6\ 2 \ cB 

The particular solutions for uniform loading are given by 

/ A 2 

*P 1 

(p2=B p x i 2 


where 


4 3>t + l x 4 -2 


X 2 +l 


1+ m . ox2 


— M 


l-M 


8 ; 1 P (2 X 2 — l)X 4 — (8X 2 — 1)(4X 2 — 1 ) 

1 — u 


p / 4 

I)c 


4 M X 2 /ox2 


l 


(2 X 2 — 1)+ 1 +4X 2 


B v = — — — 

4 8 ! (2X 2 — 1 ) X 4 — (8 X 2 — 1 )(4 X* — 1) 


p/ 4 

l)c 2 


1- 


M 


The general solution is the sum of the homogeneous solu- 
tions and the particular integral 

<p\' = ^2 A n Xi y n _I + A p X\ 

n = 1 


<P‘2 — ^2 H n Zi 7h l -\-B p x i 2 (B 1 7) 

71 = 1 

where the values y n are the six roots of the characteristic 
equation (BIG) and the coefficients A n and B n are the co- 
efficients corresponding to each of these roots. After inte- 
gration (pi becomes 

(B18) 

71 = 1 Yn o 

The general solution for <p 0 from equation (Bll) is found 
to be 

<Po= 2 2 C n Xi yn+1 -\-CpXi i -]-C <l Xi-\-Cr (B 1 9) 

71 = 1 


where, for 71 = 1 , 2, . . . (>, 


Cn = 


7n(Yn+ lj 


[ y " 2 1 + ^ [(7. - 1 )(7. - 2) + } 


and 


c Af A + 2t/ 2,X;\ pf I 

' 1 2 L ' 3 V + l-n/ ’ 6DcJ 


The coefficients Ai to A 6 , A qj C q , and (7 r must be determined 
by the boundary conditions (Bl2) to (B15). 

A complete set of coefficients is given in the following table 
for delta plates with Poisson’s ratio p. equal to 1/3 and with 


l 3 

\=-=l and Deflection curves plotted from these 

C o 

results are shown in figures 2 and 3 in which the 15° plate 


/ l 31 

corresponds to — = 1 and the 60° plate corresponds to — =\— 

C Co 


m 

7m 

4 I)C 

Am —r. 

pi* 

Pi 4 

c -5. 
m pi* 

x = 1 

x=V? 

3 

X = 1 

x=VI 

3 

x = 1 

x=V! 

3 

X = 1 

x=VI 

3 

1 

2 

3 

4 

5 

6 
V 
Q 
r 

2. 7034 
4. 9437 
8. 3816 
-2. 7034 
-4. 9437 
-8. 3816 

1. 5671 
3. 6347 
4.7258 
-1.5671 
-3. 6347 
-4. 7258 

0. 7378 
.02411 
.03827 

0 

0 

0 

-.8000 

-.01557 

0. 09632 
. 3707 
-. 1766 
0 
0 
0 

-.2903 
-. 02924 

-0. 3133 
-.03039 
-.006293 
0 
0 
0 

.3500 

-0. 1022 
-.4313 
. 07379 

0 

0 

0 

.4597 

-0. 02931 
. 003074 
. 000486 

0 

0 

0 

.04167 
-.07152 
. 05668 

-0. 003223 
.01347 
.002317 

0 

0 

0 

. 004032 
-.08354 
. 06692 


Substitution of equations (Bl7), (Bl8), and (Bl9) into the 
equation 

w=<Po J rV<P\ + 2 /V 2 


gives the expression for the deflection w of the plate under the 
assumption of parabolic chordwise deflection. 
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DELTA PLATE OF SYMMETRICAL DOUBLE- WEDGE AIRFOIL SECTION WITH 
CONSTANT THICKNESS RATIO UNDER UNIFORM LOAD 

For a delta plate of symmetrical double-wedge airfoil 
section with a constant thickness ratio the thickness is a 
function of x and y and is gi\ r en by the following equations: 

(o ^U^ C ~y) 

t= 6<«.(*i-f) (—^rjZcx^ 

where f av is the average thickness. From these expressions 
for the thickness the stiffness can be found and the coeffi- 
cients in the differential equations become 


ar- 


a>= 


21Dc 
27 Dc 2 


X\ 


8 

\)I)c* 




a.,= 5 *,• 

si 77(- 4 . 

a A- g 0 *, 


_2073 Dr' . 
4480 X ' 




\)I) c 3 
80 


« 6 
X\ 


■ _ 9Dc* . 

2 80 * 1 

. _ 26 18 TV' 

22400 X ' 

pC" n 
Vn = - ** 


Solution for linear chordwise deflections. — By use of the 
coefficients given by equations (B20), equations (A22) and 


(B20) 


(Al8) for linear chordwise deflections may be solved for <pi 
and <p 0 . The steps in the solution are the same in form as 
those for the uniform-thickness plate and the resulting 
equations are 


*>i = 


40 


3 

1 —Xi r ~2 


'^-1)1 


o +l0g,Z: § 


pi* 


7—< 


Dc 


(B21) 


and 


‘ f>0 ~27 


20 7 9 . 


y ~4 y ~ 


Xi 7 2—1 


y ~o 


+ 1— Xi 1 + 


20 2 \ 

g + 3 1 (*1 logs *1 + 1 —Xi) 


K T ~4 


C (B22) 

D 


where 




+ 80\- 


Solution for parabolic chordwise deflections. — By use of 

the coefficients given by equations (B20), equations (A41), 
(A3 5), and (A.33) for parabolic chordwise deflections may be 
solved for <p \ . <P 2 , and >po. The steps in the solution are 
again the same in form as those for the uniform-thickness 
plate and the resulting general expressions for <p h <p 2 , and 
<Po are 




j X\ 7 2 l 

— — \-A p \og e Xi 

7i = 1 O 

7«— 7T 


^2=2 B n x y " 2 +B„ — 

71 = 1 Xl 

6 1 

V 3 0=Xj C n Xi "‘ 2-| -C p X\ loge X\-\-C t Xi-\-C r 


(B23) 

(B24) 

(B25) 


where the exponents y„, are the roots of the characteristic equation 


For n— 1, 2, ... 6, A n , B n , and G„ are related by 

B n =- 


[^( 72 ~t) _ 200x2 ] =° (B26 ) 


7 2 — — 80X 2 
4 


A, 


( 7- i)[( 7 “i)( 7+ i) _80x2 ] * 


c n = 


M 

„ A 


3 ) 2 

u 2 ( 


»+T7 C 
1 D 


[( 7 --f) ( 7 "-y+;^> 
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For uniform load the coefficients in the particular integrals of equations (B23), (B24), and (B25) are 


. 10 vl* 

A P=07^ 


40 pP 

n p 9 ~ “ 


^i+1024X 2 + 320 ' 1 M X 4 -16 (5+. X 2 VlOX 2 +l) 
i 1— M ' 1— M / 

7 Dc (20X 2 +1) ^ S ‘. 1 -+ 1024X 2 +320 [ X 4 )- 1 20 (1 6X 2 + 1) (1 0X 2 + 1) 

15(16X 2 +1)-2(20X 2 +1)(5+ 1 2 ^') 

2 ‘ Dc * (20X 2 +1) (^p + 1024X 2 + 320 J ^ X 4 ^— 1 20 (1 6X 2 + 1) (1 0X 2 + 1) 

r _c , 8 2 / ,5 fi)i ! \ ,, . 2 pl^ 

C "~2 Ar 15 ( + 2I- M ) " + 81 U 

The coefficients A x to A 6 , A Qf C Q , and C r are again determined by the boundary conditions (A26), (A36), (A38), and (A39) in 
which the coefficients given by equations (B 20 ) are substituted. 

For Poisson’s ratio u equal to \ and X= -= 1 , the solution of the characteristic equation (B26) leads to two real values and 

3 c 

two pairs of complex conjugate values for 7 . The identity 

Xi a±ib—X i a cos (6 log t> X l )±ix i a sill (6 log* Xi) 

was therefore used to transform the terms involving the complex conjugate values into real form. II -= 1 and u=-> the 
solution is 

*> 2 =ip I” 0.00407Ctei 3 - 947 — 0.004363xi 8 075 cos (2.825 log, Xi)+0.006893x! 8 075 sin (2.825 log, Xi) +0.000294 


v l 


=]£- f_ o.003896x, 4 947 + 0.002 1 342+ 075 cos (2.825 log, Xi) — 0.006381x, 9 075 sin (2.825 log, xO + 0.01794 log, Xi+0.0C 1768 

Dc L -I 


<Po=?j= ^0 .00077 1 5x t 5 947 — 0 .00007 08xi 10 075 cos (2.825 log,x,) + 0.001234X, 10075 sin (2.825 log, x,) + 


0.0333 lx! log, x!-0.04096x,+ 0.040: 


26 j 


APPENDIX C 

NUMERICAL PROCEDURE FOR SOLVING DIFFERENTIAL EQUATIONS 


In cases where the equations of the present theory cannot 
be solved exactly a numerical method must be used. In 
this appendix, equations (A19) and equations (A34) and 
(A35) are set up in difference form for numerical solution. 
Initially the assumption is made that the functions involved 
in the differential equations are continuous and nonsingular. 
In this case, first and second derivatives can be expressed by 
the standard difference forms 

(d 2 y\ _y n +i—2y T ,+y n _ l 
\dx 2 )n 

(dy\ = y n +\~ yn -{ 

\dx ) n e 

where e is the distance between equally spaced station points. 

In the following development five equally spaced span- 
wise stations are used; however, the extension to a different 
number of stations can be readily made. 


The quantity (6 1 <^ l ") / at the half-stations can be expressed in 


First, consider equation (Al9) resulting from the linear 
theory 

(bm'T- 2(1 - M ) (a m 'Y=p 2 -(^£ J‘ Pl (x) dx 2 J -g, 

Because of the nature of the tip boundary conditions for 
this equation, it can be conveniently put in the form 

T' = q } (Cl) 

where 

T= (b]#/')' —2(1— n)a } (pi' 


In finding the difference equation equivalent to equation 
(Cl), the quantity (bupi")' is found in matrix form; from 
this expression is subtracted the matrix equivalent of 
2(1— }i)ai<pi) the resulting expression for T is multiplied by 
a differentiating matrix; and the product is equated to the 
right-hand side. 

matrix form as follows: 


rr 

V 7 ! 0 


~i 

-2 1 


“ 

rr 

<Pl 1 



1 -2 1 



rr 

1 2 

_ 1 

2 


1 -2 

1 



e 





rr 

&\ 3 



1 

-2 1 


rr 

<Pl 4 




1 -2 

1_ 


<pi-i 

<Pio 

<Pn 

*P\2 

<Pl3 

<PU 

^15 


(C2) 


where the second subscript denotes the station point, the subscript at the root station being 0 and at the tip 5. The 
root boundary conditions are now applied; namely, 


<Pi ( 0 ) = 0=^,0 


r (f\\ ft l ^l—l 

<p i (0)=o= — Y e — 

Thus, after the values of <p\o—f) and <pi-i = <pn are substituted, equation (C2) becomes 


<P\'q 


”2 


<Pu 

<P\"\ 

1 

e 2 

-2 1 


<Pl2 

rr 

<P 1 2 

rr 

<P 1 3 

1 -2 1 
1 -2 1 


<Pl3 

<P 14 

rr 

<P 1 4 


1 —2 1_ 


<Pl5 


15 


1G 

Therefore, 
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(6,¥>i")o 


(6,<p,")i 


(WO* 

1 

~ 7 2 

(6m*). 





6,1 


b V 2 


b\3 


bu. 


2 

-2 1 
1 -2 
1 


1 

-2 1 
1 -2 


<Pu 

<Pl2 

<Pl3 

<£>14 

<Pl5 


One of the tip boundary conditions is 
Thus, 


(bi<pi") X m i — 0 — (bi(pi") 5 


//^/ 

1/2 


(b m") 
(bm'%, 
(bm”) 
(b\(p\") 7 /2 

//\f 


/t\f 

5/2 


— 1 1 
-1 


(bi(p\") Qf 2 

The matrix equivalent for the second term of T is 


1 

1 1 
-1 1 
- 1 . 


(bm") 

(&i*0 

(bi<pi") 

{b m *) A 


2(1— /x) 


(d\(fl')l/2 
(dl<Pi f )s/2 
(di<Pi f ) 5/2 
(ai^i') 7 / 2 
(d\<P\')Q/2 


2(1 — m) 


*- 1 , 1/2 


*1,3/2 


d], 5/2 


d},7/2 


*'1,9/2 —I 


‘ 1 

-1 1 
— 1 


1 

— 1 


Therefore T becomes 


Tl/2 

T 3/2 

^5/2 

Ti/2 

T9/2 


1 

-1 1 
— 1 


-1 


2(1— M ) 


&i, 1/2 


* 1 , 3/2 


1 

-1 


*-1,5/2 


&11 


6l2 


*1,7/2 


&1, 9/2- 


b\s 

b 

" 1 

— 1 


2 

— 2 1 
1 -2 1 
1 -2 


1 -2 


1 

-1 


1 

-1 


<Pn 

^12 

<Pl3 

<Pl4 

<Pl5 


(C3) 


(C4) 


1 1 


<Pn 

<P 12 

<Pl3 

^14 

¥>15 


The right-hand side of equation (Cl) can now be equated to the derivative of equation (C5) ; thus, 


2n 

</l2 

#13 

#14 


-1 1 
— 1 


1 

— 1 


T 1/2 
Tz/2 
Tz/2 

T7/2 

T9/2 


In order to obtain gi 5 , the boundary condition 


T= 0 


(C5) 


at x 


c=i must be used. In other words, T goes from T m at station to 0 at station 5. A straight line drawn between 
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9 J 7 

these two points would have the slope — - J' 2 - The value of q 15 is considered to be this slope; therefore.. 


$11 


"—11 


T \/2 

$12 

1 

-1 1 


Tz /2 

$13 

€ 

-1 1 


T $/2 

$14 


-1 1 


Tt /2 

2 . 5/2 


— 1. 


T 9 /2 


or 


$11 



"1-2 1 

610 

”2 

$12 


1 

1 -2 1 

bn 

-2 1 

$13 

= < 

i 

? 

1 -2 1 

b \ 2 

1 -2 1 

$14 



1 -2 

b\z 

1 -2 1 

9.5/2 



1. 

bu_ 

1 -2 1_ 


1 1 

#i, 1/2 

" 1 

-1 1 

# 1 , 3/2 

-1 1 

-1 1 

# 1 , 5/2 

— 1 1 

-1 1 

# 1 , 7/2 

-1 1 

— 1_ 

# 1 , 9 / 2 . 

— 1 


If the matrix multiplication is carried out, the difference equivalent of equation (Cl) finally becomes 


$11 


<Pll 

$12 


<Pl 2 

$13 


<£>13 

$14 


<£>14 

9.5/2 


<Plo 


where 



26io + 46ii + 6i2 

— 26n — 26i2 

612 


- 


26n — 26i2 

6ll + 46i2 + 6i3 

26i2 — 26i3 

613 


[C,] = 


612 

26i2 26i3 

612+4613 + 6 14 

2613 2614 

bu 



613 

2613 26 h 

613+4614 

26l4 




614 

-2614 

6l4 _ 


” . 

#1, 1/2 — #1,3/2 

#1,3/2 






#1,3/2 — 

#1,3/2 — #1,5/2 

#1,5/2 



[£>.1= 



#1,5/2 

#1,5/2 #1,7/2 

#1,7/2 






#1,7/2 

#1,7/2 “#1,9/2 

#1,9/2 


_ 




#1,9/2 

#1,9/2 


In order to determine <p 0 from <p u use must be made of equation (Al8) 



Pi dx 2 —— <p\ 

# i 


n 


or, by use of the boundary condition <p 0 (0) = «/(0)=0, 


<Po — 






Pi 


dx A — 



a 2 


■ <p i 


dx 2 


IJJ 


<Pn 

<PV 2 

^13 

<Pl 4 

<Pl 5 


(C6) 


(C7) 
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In matrix form equation (C7) becomes 


<A)i 

<£>02 

<^03 

<A>4 


1 1 
1 1 
1 1 
1 1 


1 1 


1 1 1 
1111 


1/2 
1/2 
1/2 
1/2 

Ll/2 

■ 1/2 
1/2 1 
1/2 1 1 
1/2111 


1/011 


1/a, 


1/0,2 


1/013 


I/01J 


1 1 1 lj L 1/2 1 1 1 lj _ 


'(02/01)0 


(02/01)1 


(02/0O2 


(02/01)3 


'1111 
1 1 1 
1 1 


9 


1 -2 


1 -2 


P11 

P12 

Pl3 

Pu 

Pu/1 

<P \‘2 
<P\Z 
<Pu 


(a 2 /ai)J _ 


1 -2 lj <p l5 


(C8) 


Thus, if the values of q\ (which can be determined numerically or analytically according to preference and feasibility) are 
known, the values of <p v can be found by solving equation (C6) and the values of <p 0 in turn by means of equation (C8). 

The foregoing development applies to the case where only linear chordwise deformations are allowed. A similar pro- 
cedure is followed in expressing the differential equations pertaining to the parabolic theory in difference form; only the 
results are shown herein. 


The matrix equivalent to equations (A34) and (A35) is 


<lu 

q \2 

</l3 

<Zl4 

<Zio/2 


p [C’.]- 2(1 7 — [DA 

7 [DA 


<£ll 

<Pl 2 

<Pu 

<Pl 5 

<221 





<P 2 l 

q22 





<P 22 

fe 


-4 [DA 

6 € 

jACA-^t^[DA+ 


<£>23 

tf24 



4(1-m 2 )[A’1 


<P 24 

<1 25/2 





<P2o 


where 



2b n0 -{-4:b n i-\-b n2 

—2b nl —2b n2 

b n 2 





— 2b nX —2b n2 

^nlH“46 n 2+6»3 

— 2b n2 26 /;3 

b n Z 



[Cnl- 


b n 2 

—2b n2 —2b n 3 

ft«2+46 n3 +6„4 

— 26 /i 3 — 26 w4 

b nA 




bnZ 

— 2b n3 —2b ni 

&ti3 + 46„4 

~2b n < 





bnA 

-26*4 

bnA — 



- 

“ #w, 1/2 — #«,3/2 

# 71 , 3/2 




' 



# 77 , 3/2 

— #71,3/2 # 77 , 5/2 

#w,5/2 



m= 



#/i, 5/2 

— #71,5/2 # w, 7/2 

#71,7/2 






#7»,7/2 

— #71,7/2 — #71,9/2 

#71,9/2 






# 77 , 9/2 

#71,9/2. 


(C9) 
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®ll 


a 12 


[£] = 


and qi and q 2 are the right-hand sides of equations (A34) and (A35), respectively; that is, 

? ' =p, “(S ],'/>■ -)' 

q ,-p,-2 „££„ ££p, dx ■)' 

With and <p 2 known, <po can be obtained by use of equation (A33) 

1 / 10 


^01 

<P02 

<A>3 

<A)4 

^05 


1 

1 1 
1 1 1 
1111 
11111 


' 1/2 
1/2 1 
1/2 1 1 
1/2111 
. 1/2 1111 . 


l/dn 


l/a l2 


1/tf 13 


1/a i4 J 


11111“ 
1111 
1 1 1 
1 1 
1 


11111" 
1111 
1 1 1 
1 1 
1 


Vn 

Pl2 

Vn 

Pu 

Pn/2 


‘(a 2 /ai) 0 


(a 2 /a l ) l 


(a 2 /ai) 2 


(a 2 /ai) 3 


fe/ai)4- 


(a 3 /ai) 0 


(a 3 /ai)i 


(a 3 /ai) 2 


(a 3 /ai) 3 


(a 3 /ai) 4 _ 


■ 2 

-2 1 

1 —2 1 

1 -2 1 

1 -2 1. 

■ 2 

-2 1 

1 —2 1 

1 -2 1 

1 -2 1 


<Pn 

<Pl2 

<Pl3 

<Pu 

<Pl5 

<P‘2l 

<P-22 

*P'2Z 

<£>24 

^25 




<?21 

<£>22 

<P23 

<?24 

^25 


(CIO) 


It should be noted that, as can be expected, the matrix 
equations (C6) and (C8) are merely special cases of equations 
(C9) and (CIO), respectively. In addition, the square 
matrices in equations (C6) and (C9) are symmetric, a result 
that is consistent with the fact that the differential equations 
under consideration are self-adjoint. 

In the beginning of this appendix the assumption was made 
that the functions involved in the differential equations are 
continuous and nonsingular. The difference solution, how- 
ever, may be adequate for some cases in which this assump- 
tion is not strictly correct. For instance, the deflections of 
a plate with a discontinuous stiffness distribution coidd con- 
ceivably be not very different from the deflections of a plate 


with a continuous stillness distribution closely approximating 
the discontinuous distribution except in the neighborhood of 
the discontinuity. The results yielded by the difference solu- 
tion in this case would be those associated with the continuous 
stiffness distribution. The number of stations may have to 
be increased, however, in order to minimize the inaccuracy 
introduced by the discontinuity or, in other cases, by a 
singularity. The case of the symmetrical double-wedge air- 
foil section, constant-thickness-ratio delta plate, discussed 
in the body of this report, is an example of a treatment of a 
singularity. In this case, although the solution is singular, 
adequate accuracy is obtained by the difference solution if 
ten equal intervals are used. 


APPENDIX D 


DEFLECTION AND STRESS EXPERIMENTS 

Test specimens. — The specimens tested were: (1) a 45° 
right-triangular plate clamped along one leg and (2) a 60° 
right-triangular plate clamped along the longer leg. Each 
specimen, cut from 24S-T4 aluminum-alloy sheet of 0.250- 
inch thickness, had a length perpendicular to the clamped 
edge of 30 inches. 

Method of testing. -Figure 12, a photograph of the test 
setup, shows the methods of clamping, loading, and measure- 
ment of deflections. A 1,000,000-pound clamping load (held 
constant during the test) was applied to the root area of each 
specimen and a uniform load of 0.204 psi was applied by 
2-incli washers giving a tip deflection in each case of approx- 
imately £ inch. 

The deflections were measured by dial gages placed at the 
points indicated in figures 2 and 3. 

Stresses were obtained from the 45° specimen only. On 
this specimen, 13 resistance-wire rosette strain gages were 
placed at the points indicated in figure 5. The plate was 
loaded with 2-inch washers in four increments of 0.0847 psi 
per increment and the maximum tip deflection was 1.13 
inches. Readings of all the strain gages were recorded at 
each increment of loading. 

Analysis and discussion of data. — The deflection w was 
plotted in figures 2 and 3 in terms of the nondimensional 



Figure 12. — Deflection test setup of the 45° delta plate under uniform 

load. 
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parameter wD/pl *, in which the elastic constants were taken 
as i?=10.6X10 6 psi and It was found that the dial- 

o 

gage forces reduced the tip deflection of the plate by approx- 
imately 2 percent ; however, since this error is of ihe same 
order of magnitude as that in the material properties and 
from other sources, no corrections are made in the results 
presented. 

The readings of each of the 39 individual strain gages were 
plotted against load, and the slope of each of the resulting 
linear curves was taken as the average strain per unit load 
of the individual gage. The principal stresses were then 
calculated and plotted in figure 5 in terms of the nondimen- 
sional parameter at 2 lpl 2 . 
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